Range-wide genetic analysis of an endangered bumble bee (Bombus affinis, Hymenoptera: Apidae) reveals population structure, isolation by distance, and low colony abundance.

Declines in bumble bee species range and abundances are documented across multiple continents and have prompted the need for research to aid species recovery and conservation. The rusty patched bumble bee (Bombus affinis) is the first federally listed bumble bee species in North America. We conducted a range-wide population genetics study of B. affinis from across all extant conservation units to inform conservation efforts. To understand the species' vulnerability and help establish recovery targets, we examined population structure, patterns of genetic diversity, and population differentiation. Additionally, we conducted a site-level analysis of colony abundance to inform prioritizing areas for conservation, translocation, and other recovery actions. We find substantial evidence of population structuring along an east-to-west gradient. Putative populations show evidence of isolation by distance, high inbreeding coefficients, and a range-wide male diploidy rate of ~15%. Our results suggest the Appalachians represent a genetically distinct cluster with high levels of private alleles and substantial differentiation from the rest of the extant range. Site-level analyses suggest low colony abundance estimates for B. affinis compared to similar datasets of stable, co-occurring species. These results lend genetic support to trends from observational studies, suggesting that B. affinis has undergone a recent decline and exhibit substantial spatial structure. The low colony abundances observed here suggest caution in overinterpreting the stability of populations even where B. affinis is reliably detected interannually. These results help delineate informed management units, provide context for the potential risks of translocation programs, and help set clear recovery targets for this and other threatened bumble bee species.

Metabolomes of bumble bees reared in common garden conditions suggest constitutive differences in energy and toxin metabolism across populations.

Cold tolerance of ectotherms can vary strikingly among species and populations. Variation in cold tolerance can reflect differences in genomes and transcriptomes that confer cellular-level protection from cold; additionally, shifts in protein function and abundance can be altered by other cellular constituents as cold-exposed insects often have shifts in their metabolomes. Even without a cold challenge, insects from different populations may vary in cellular composition that could alter cold tolerance, but investigations of constitutive differences in metabolomes across wild populations remain rare. To address this gap, we reared Bombus vosnesenskii queens collected from Oregon and California (USA) that differ in cold tolerance (CTmin = -6 °C and 0 °C, respectively) in common garden conditions, and measured offspring metabolomes using untargeted LC-MS/MS. Oregon bees had higher levels of metabolites associated with carbohydrate (sorbitol, lactitol, maltitol, and sorbitol-6-phosphate) and amino acid (hydroxyproline, ornithine, and histamine) metabolism. Exogenous metabolites, likely derived from the diet, also varied between Oregon and California bees, suggesting population-level differences in toxin metabolism. Overall, our results reveal constitutive differences in metabolomes for bumble bees reared in common garden conditions from queens collected in different locations despite no previous cold exposure.

It is buzziness time: rearing, mating, and overwintering Bombus vosnesenskii (Hymenoptera: Apidae).

Bombus vosnesenskii Radowszkowski, 1862 is one of three bumble bee species commercially available for pollination services in North America; however, little is documented about B. vosnesenskii colony life cycle or the establishment of ex situ rearing, mating, and overwintering practices. In this study, we documented nest success, colony size, and gyne production; recorded the duration of mating events; assessed overwintering survival of mated gynes; and evaluated second-generation nest success for colonies established from low- and high-elevation wild-caught B. vosnesenskii gynes. Of the 125 gynes installed, 62.4% produced brood cells (nest initiation) and 43.2% had at least 1 worker eclose (nest establishment). High-elevation B. vosnesenskii gynes had significantly higher nest initiation and establishment success than low-elevation gynes. However, low-elevation colonies were significantly larger with queens producing more gynes on average. Mating was recorded for 200 low-elevation and 37 high-elevation gynes, resulting in a mean duration of 62 and 51 min, respectively. Mated gynes were then placed into cold storage for 54 days to simulate overwintering, which resulted in 59.1% of low-elevation gynes surviving and 91.9% of high-elevation gynes surviving. For second-generation low-elevation gynes, 26.4% initiated nesting and 14.3% established nesting. Second-generation high-elevation gynes did not initiate nesting despite CO2 narcosis treatments. Overall, these results increase our understanding of B. vosnesenskii nesting, mating, and overwintering biology from 2 elevations. Furthermore, this study provides information on successful husbandry practices that can be used by researchers and conservationists to address knowledge gaps and enhance the captive rearing of bumble bees.

Variation in North American bumble bee nest success and colony sizes under captive rearing conditions.

Of the 265 known bumble bee (Bombus) species, knowledge of colony lifecycle is derived from relatively few species. As interest in Bombus commercialization and conservation grows, it is becoming increasingly important to understand colony growth dynamics across a variety of species since variation exists in nest success, colony growth, and reproductive output. In this study, we reported successful nest initiation and establishment rates of colonies and generated a timeline of colony development for 15 western North American Bombus species, which were captively reared from wild-caught gynes from 2009 to 2019. Additionally, we assessed variation in colony size among 5 western North American Bombus species from 2015 to 2018. Nest initiation and establishment rates varied greatly among species, ranging from 5-76.1% and 0-54.6%, respectively. Bombus griseocollis had the highest rates of nest success across the 11-yr period, followed by B. occidentalis, B. vosnesenskii, and B. huntii. Furthermore, days to nest initiation and days to nest establishment varied among species, ranging from 8.4 to 27.7 days and 32.7 to 47 days. Colony size also differed significantly among species with B. huntii and B. vosnesenskii producing more worker/drone cells than B. griseocollis, B. occidentalis, and B. vancouverensis. Additionally, gyne production differed significantly among species with B. huntii colonies producing more gynes than B. vosnesenskii. Results from this study increase knowledge of systematic nesting biology for numerous western North American Bombus species under captive rearing conditions, which can further improve rearing techniques available to conservationists and researchers.

Comparative analysis of 3 pollen sterilization methods for feeding bumble bees.

Pollen is an essential component of bee diets, and rearing bumble bees (Bombus spp.) for commercial use necessitates feeding pollen in mass quantities. This pollen is collected from honey bee (Apis mellifera L.) colonies because neither an artificial diet nor an economical, large-scale pollen collection process from flowers is available. The provenance of honey bee-collected pollen is often unknown, and in some cases has crossed international borders. Both deformed wing virus (DWV) and the fungal pathogen Ascosphaera apis (Claussen) Olive & Spiltoir (cause of chalkbrood disease); occur in honey bee-collected pollen, and infections have been observed in bumble bees. We used these pathogens as general surrogates for viruses and spore-forming fungal diseases to test the efficacy of 3 sterilization methods, and assessed whether treatment altered pollen quality for the bumble bee. Using honey bee-collected pollen spiked with known doses of DWV and A. apis, we compared gamma irradiation (GI), ozone fumigation (OZ), and ethylene oxide fumigation (EO) against an untreated positive control and a negative control. Following sterilization treatments, we tested A. apis spore viability, detected viral presence with PCR, and tested palatability to the bumble bee Bombus impatiens Cresson. We also measured bacterial growth from pollens treated with EO and GI. GI and EO outperformed OZ treatment in pathogen suppression. EO had the highest sterilizing properties under commercial conditions and retained palatability and supported bee development better than other treatments. These results suggest that EO sterilization reduces pathogen risks while retaining pollen quality as a food source for rearing bumble bees.

Whole genome demographic models indicate divergent effective population size histories shape contemporary genetic diversity gradients in a montane bumble bee.

Understanding historical range shifts and population size variation provides an important context for interpreting contemporary genetic diversity. Methods to predict changes in species distributions and model changes in effective population size (N e) using whole genomes make it feasible to examine how temporal dynamics influence diversity across populations. We investigate N e variation and climate-associated range shifts to examine the origins of a previously observed latitudinal heterozygosity gradient in the bumble bee Bombus vancouverensis Cresson (Hymenoptera: Apidae: Bombus Latreille) in western North America. We analyze whole genomes from a latitude-elevation cline using sequentially Markovian coalescent models of N e through time to test whether relatively low diversity in southern high-elevation populations is a result of long-term differences in N e. We use Maxent models of the species range over the last 130,000 years to evaluate range shifts and stability. N e fluctuates with climate across populations, but more genetically diverse northern populations have maintained greater N e over the late Pleistocene and experienced larger expansions with climatically favorable time periods. Northern populations also experienced larger bottlenecks during the last glacial period, which matched the loss of range area near these sites; however, bottlenecks were not sufficient to erode diversity maintained during periods of large N e. A genome sampled from an island population indicated a severe postglacial bottleneck, indicating that large recent postglacial declines are detectable if they have occurred. Genetic diversity was not related to niche stability or glacial-period bottleneck size. Instead, spatial expansions and increased connectivity during favorable climates likely maintain diversity in the north while restriction to high elevations maintains relatively low diversity despite greater stability in southern regions. Results suggest genetic diversity gradients reflect long-term differences in N e dynamics and also emphasize the unique effects of isolation on insular habitats for bumble bees. Patterns are discussed in the context of conservation under climate change.

Recent and future declines of a historically widespread pollinator linked to climate, land cover, and pesticides.

The acute decline in global biodiversity includes not only the loss of rare species, but also the rapid collapse of common species across many different taxa. The loss of pollinating insects is of particular concern because of the ecological and economic values these species provide. The western bumble bee (Bombus occidentalis) was once common in western North America, but this species has become increasingly rare through much of its range. To understand potential mechanisms driving these declines, we used Bayesian occupancy models to investigate the effects of climate and land cover from 1998 to 2020, pesticide use from 2008 to 2014, and projected expected occupancy under three future scenarios. Using 14,457 surveys across 2.8 million km2 in the western United States, we found strong negative relationships between increasing temperature and drought on occupancy and identified neonicotinoids as the pesticides of greatest negative influence across our study region. The mean predicted occupancy declined by 57% from 1998 to 2020, ranging from 15 to 83% declines across 16 ecoregions. Even under the most optimistic scenario, we found continued declines in nearly half of the ecoregions by the 2050s and mean declines of 93% under the most severe scenario across all ecoregions. This assessment underscores the tenuous future of B. occidentalis and demonstrates the scale of stressors likely contributing to rapid loss of related pollinator species throughout the globe. Scaled-up, international species-monitoring schemes and improved integration of data from formal surveys and community science will substantively improve the understanding of stressors and bumble bee population trends.

Novel Microsatellite Markers for Osmia lignaria (Hymenoptera: Megachilidae): A North American Pollinator of Agricultural Crops and Wildland Plants.

Comprehensive decisions on the management of commercially produced bees, depend largely on associated knowledge of genetic diversity. In this study, we present novel microsatellite markers to support the breeding, management, and conservation of the blue orchard bee, Osmia lignaria Say (Hymenoptera: Megachilidae). Native to North America, O. lignaria has been trapped from wildlands and propagated on-crop and used to pollinate certain fruit, nut, and berry crops. Harnessing the O. lignaria genome assembly, we identified 59,632 candidate microsatellite loci in silico, of which 22 were tested using molecular techniques. Of the 22 loci, 12 loci were in Hardy-Weinberg equilibrium (HWE), demonstrated no linkage disequilibrium (LD), and achieved low genotyping error in two Intermountain North American wild populations in Idaho and Utah, USA. We found no difference in population genetic diversity between the two populations, but there was evidence for low but significant population differentiation. Also, to determine if these markers amplify in other Osmia, we assessed 23 species across the clades apicata, bicornis, emarginata, and ribifloris. Nine loci amplified in three species/subspecies of apicata, 22 loci amplified in 11 species/subspecies of bicornis, 11 loci amplified in seven species/subspecies of emarginata, and 22 loci amplified in two species/subspecies of ribifloris. Further testing is necessary to determine the capacity of these microsatellite loci to characterize genetic diversity and structure under the assumption of HWE and LD for species beyond O. lignaria. These markers will inform the conservation and commercial use of trapped and managed O. lignaria and other Osmia species for both agricultural and nonagricultural systems.

Captive Rearing Success and Critical Thermal Maxima of Bombus griseocollis (Hymenoptera: Apidae): A Candidate for Commercialization?

Commercialized bumble bees (Bombus) are primary pollinators of several crops within open field and greenhouse settings. However, the common eastern bumble bee (Bombus impatiens Cresson, 1863) is the only species widely available for purchase in North America. As an eastern species, concerns have been expressed over their transportation outside of their native range. Therefore, there is a need to identify regionally appropriate candidates for commercial crop pollination services, especially in the western U.S.A. In this study, we evaluated the commercialization potential of brown-belted bumble bees (Bombus griseocollis De Geer, 1773), a broadly distributed species throughout the U.S.A., by assessing nest initiation and establishment rates of colonies produced from wild-caught gynes, creating a timeline of colony development, and identifying lab-reared workers' critical thermal maxima (CTMax) and lethal temperature (ecological death). From 2019 to 2021, 70.6% of the wild-caught B. griseocollis gynes produced brood in a laboratory setting. Of these successfully initiated nests, 74.8% successfully established a nest (produced a worker), providing guidance for future rearing efforts. Additionally, lab-reared workers produced from wild-caught B. griseocollis gynes had an average CTMax of 43.5°C and an average lethal temperature of 46.4°C, suggesting B. griseocollis can withstand temperatures well above those commonly found in open field and greenhouse settings. Overall, B. griseocollis should continue to be evaluated for commercial purposes throughout the U.S.A.

Bumble bees exhibit body size clines across an urban gradient despite low genetic differentiation.

Environmental heterogeneity resulting from human-modified landscapes can increase intraspecific trait variation. However, less known is whether such phenotypic variation is driven by plastic or adaptive responses to local environments. Here, we study five bumble bee (Apidae: Bombus) species across an urban gradient in the greater Saint Louis, Missouri region in the North American Midwest and ask: (1) Can urban environments induce intraspecific spatial structuring of body size, an ecologically consequential functional trait? And, if so, (2) is this body size structure the result of plasticity or adaptation? We additionally estimate genetic diversity, inbreeding, and colony density of these species-three factors that affect extinction risk. Using ≥ 10 polymorphic microsatellite loci per species and measurements of body size, we find that two of these species (Bombus impatiens, Bombus pensylvanicus) exhibit body size clines across the urban gradient, despite a lack of population genetic structure. We also reaffirm reports of low genetic diversity in B. pensylvanicus and find evidence that Bombus griseocollis, a species thought to be thriving in North America, is inbred in the greater Saint Louis region. Collectively, our results have implications for conservation in urban environments and suggest that plasticity can cause phenotypic clines across human-modified landscapes.

Neonicotinoid Pesticides Cause Mass Fatalities of Native Bumble Bees: A Case Study From Wilsonville, Oregon, United States.

In June of 2013 an application of dinotefuran on an ornamental planting of European linden trees (Tilia cordata Mill. [Malvales: Malvalceae]) in a shopping mall parking lot in Wilsonville, Oregon provoked the largest documented pesticide kill of bumble bees in North America. Based on geographic information systems and population genetic analysis, we estimate that between 45,830 and 107,470 bumble bees originating from between 289 and 596 colonies were killed during this event. Dinotefuran is a neonicotinoid that is highly effective in exterminating and/or harming target pest insects and non-target beneficial insects. Analysis to detect the concentration of pesticides in flowers that received foliar application revealed that the minimum reported dinotefuran concentration of a sampled T. cordata flower was 7.4 ppm, or in excess of 737% above the LC50 of the beneficial pollinator, the honey bee (Apis mellifera Linnaeus, 1758 [Hymenoptera: Apidae]). Furthermore, sampled Vosnesensky bumble bees (Bombus vosnesenskii Radoskowski, 1862 [Hymenoptera: Apidae]) were found to have an average dinotefuran concentration of 0.92 ppm at the time of death, which exceeds the maximum LC50 of A. mellifera (0.884 ppm). Our study underscores the lethal impact of the neonicotinoid pesticide dinotefuran on pollinating insect populations in a suburban environment. To our knowledge, the documentation and impact of pesticide kills on wild populations of beneficial insects has not been widely reported in the scientific literature. It is likely that the vast majority of mass pesticide kills of beneficial insects across other environments go unnoticed and unreported.

Body mass and sex, not local climate, drive differences in chill coma recovery times in common garden reared bumble bees.

The time required to recover from cold exposure (chill coma recovery time) may represent an important metric of performance and has been linked to geographic distributions of diverse species. Chill coma recovery time (CCRT) has rarely been measured in bumble bees (genus Bombus) but may provide insights regarding recent changes in their distributions. We measured CCRT of Bombus vosnesenskii workers reared in common garden laboratory conditions from queens collected across altitude and latitude in the Western United States. We also compared CCRTs of male and female bumble bees because males are often overlooked in studies of bumble bee ecology and physiology and may differ in their ability to respond to cold temperatures. We found no relationship between CCRT and local climate at the queen collection sites, but CCRT varied significantly with sex and body mass. Because differences in the ability to recover from cold temperatures have been shown in wild-caught Bombus, we predict that variability in CCRT may be strongly influenced by plasticity.

Biogeographic parallels in thermal tolerance and gene expression variation under temperature stress in a widespread bumble bee.

Global temperature changes have emphasized the need to understand how species adapt to thermal stress across their ranges. Genetic mechanisms may contribute to variation in thermal tolerance, providing evidence for how organisms adapt to local environments. We determine physiological thermal limits and characterize genome-wide transcriptional changes at these limits in bumble bees using laboratory-reared Bombus vosnesenskii workers. We analyze bees reared from latitudinal (35.7-45.7°N) and altitudinal (7-2154 m) extremes of the species' range to correlate thermal tolerance and gene expression among populations from different climates. We find that critical thermal minima (CTMIN) exhibit strong associations with local minimums at the location of queen origin, while critical thermal maximum (CTMAX) was invariant among populations. Concordant patterns are apparent in gene expression data, with regional differentiation following cold exposure, and expression shifts invariant among populations under high temperatures. Furthermore, we identify several modules of co-expressed genes that tightly correlate with critical thermal limits and temperature at the region of origin. Our results reveal that local adaptation in thermal limits and gene expression may facilitate cold tolerance across a species range, whereas high temperature responses are likely constrained, both of which may have implications for climate change responses of bumble bees.

A 3D Printed Pollen Trap for Bumble Bee (Bombus) Hive Entrances.

To verify the plant sources from which bumble bees forage for pollen, individuals must be collected to remove their corbicular pollen loads for analysis. This has traditionally been done by netting foragers at nest entrances or on flowers, chilling the bees on ice, and then removing the pollen loads from the corbiculae with forceps or a brush. This method is time and labor intensive, may alter normal foraging behavior, and can result in stinging incidents for the worker performing the task. Pollen traps, such as those used on honey bee hives, collect pollen by dislodging corbicular pollen loads from the legs of workers as they pass through screens at the nest entrance. Traps can remove a large quantity of pollen from returning forager bees with minimal labor, yet to date no such trap is available for use with bumble bee colonies. Workers within a bumble bee colony can vary in size making size selection of entrances difficult to adapt this mechanism to commercially reared bumble bee hives. Using 3D printing design programs, we created a pollen trap that successfully removes the corbicular pollen loads from the legs of returning bumble bee foragers. This method significantly reduces the amount of time required by researchers to collect pollen from bumble bee foragers returning to the colony. We present the design, results of pollen removal efficiency tests, and suggest areas of modifications for investigators to adapt traps to a variety of bumble bee species or nest box designs.

De Novo Genome Assemblies for Three North American Bumble Bee Species: Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii.

Bumble bees are ecologically and economically important insect pollinators. Three abundant and widespread species in western North America, Bombus bifarius, Bombus vancouverensis, and Bombus vosnesenskii, have been the focus of substantial research relating to diverse aspects of bumble bee ecology and evolutionary biology. We present de novo genome assemblies for each of the three species using hybrid assembly of Illumina and Oxford Nanopore Technologies sequences. All three assemblies are of high quality with large N50s (> 2.2 Mb), BUSCO scores indicating > 98% complete genes, and annotations producing 13,325 - 13,687 genes, comparing favorably with other bee genomes. Analysis of synteny against the most complete bumble bee genome, Bombus terrestris, reveals a high degree of collinearity. These genomes should provide a valuable resource for addressing questions relating to functional genomics and evolutionary biology in these species.

Local adaptation across a complex bioclimatic landscape in two montane bumble bee species.

Understanding evolutionary responses to variation in temperature and precipitation across species ranges is of fundamental interest given ongoing climate change. The importance of temperature and precipitation for multiple aspects of bumble bee (Bombus) biology, combined with large geographic ranges that expose populations to diverse environmental pressures, make these insects well-suited for studying local adaptation. Here, we analyzed genome-wide sequence data from two widespread bumble bees, Bombus vosnesenskii and Bombus vancouverensis, using multiple environmental association analysis methods to investigate climate adaptation across latitude and altitude. The strongest signatures of selection were observed in B. vancouverensis, but despite unique responses between species for most loci, we detected several shared responses. Genes relating to neural and neuromuscular function and ion transport were especially evident with respect to temperature variables, while genes relating to cuticle formation, tracheal and respiratory system development, and homeostasis were associated with precipitation variables. Our data thus suggest that adaptive responses for tolerating abiotic variation are likely to be complex, but that several parallels among species can emerge even for these complex traits and landscapes. Results provide the framework for future work into mechanisms of thermal and desiccation tolerance in bumble bees and a set of genomic targets that might be monitored for future conservation efforts.

Why Are Queens Broodless? Failed Nest Initiation Not Linked to Parasites, Mating Status, or Ovary Development in Two Bumble Bee Species of Pyrobombus (Hymenoptera: Apidae: Bombus).

Bumble bees (Bombus [Hymenoptera: Apidae]) are important pollinators for agricultural crops, which has led to their commercial domestication. Despite their importance, little is known about the reproductive biology of bumble bees native to North America. The Hunt bumble bee (Bombus huntii Greene [Hymenoptera: Apidae]) and the Vosnesensky bumble bee (Bombus vosnesenskii Radoszkowski [Hymenoptera: Apidae] are native candidates for commercial production in western North America due to their efficacy in providing commercial pollination services. Availability of pollinators native to the region in which services would be provided would minimize the likelihood of introducing exotic species and spreading novel disease. Some parasites are known to affect bumble bee reproduction, but little is known about their prevalence in North America or how they affect queen success. Only 38% of wild-caught B. huntii and 51% wild-caught B. vosnesenskii queens collected between 2015 and 2017 initiated nests in the laboratory. Our objective was to identify causal factors leading to a queen's inability to oviposit. To address this, we dissected each broodless queen and diagnosed diseases, assessed mating status, and characterized ovary development. Nematodes, arthropods, and microorganisms were detected in both species. Overall, 20% of queens were infected by parasites, with higher rates in B. vosnesenskii. Over 95% of both species were mated, and over 88% had developed ovaries. This suggests that parasitism and mating status were not primary causes of broodlessness. Although some failure to nest can be attributed to assessed factors, additional research is needed to fully understand the challenges presented by captive rearing.

Bombus (Hymenoptera: Apidae) Microcolonies as a Tool for Biological Understanding and Pesticide Risk Assessment.

Bumble bees provide valuable pollination services to many wild and agricultural plants. Populations of some bumble bee species are in decline, prompting the need to better understand bumble bee biology and to develop methodologies for assessing the effects of environmental stressors on these bees. Use of bumble bee microcolonies as an experimental tool is steadily increasing. This review closely examines the microcolony model using peer-reviewed published literature identified by searching three databases through November 2018. Microcolonies have been successfully used for investigating a range of endpoints including behavior, the gut microbiome, nutrition, development, pathogens, chemical biology, and pesticides/xenobiotics. Methods for the initiation and monitoring of microcolonies, as well as the recorded variables were catalogued and described. From this information, we identified a series of recommendations for standardizing core elements of microcolony studies. Standardization is critical to establishing the foundation needed to support use of this model for biological response investigations and particularly for supporting use in pesticide risk assessment.

Foraging Economics of the Hunt Bumble Bee, a Viable Pollinator for Commercial Agriculture.

Globally, there are only five bumble bee (Hymenoptera: Apidae, Bombus) species that have been successfully commercialized for agriculture. The Hunt bumble bee, Bombus huntii Green, 1860, has been recognized as a suitable pollinator of crops and has a broad distribution in western North America, making it a viable candidate for commercialization. In this study, our goal was to characterize the foraging dynamics of B. huntii female workers under open field conditions. To accomplish this goal, we monitored three B. huntii colonies over an 8-wk period in the summer of 2012 in northern Utah. Using marked bees, we studied the relationship between foraging duration/offloading and pollen/nonvisible pollen collection. In total, we observed 921 foraging events across all three colonies. Of our observations, 82% (n = 756) were foraging events that included both a departure and arrival time observation. Average duration of pollen and nonpollen (i.e., nectar) trips across foragers is 41.86 ± 5.65 min (±SE) and 32.18 ± 5.89 min, respectively. Workers spent a significantly longer time offloading pollen in the nest after a foraging trip relative to workers without pollen present on their corbicula. Pollen foraging rate increases over the course of the day, likely due to the time it takes to learn how to forage on a diverse array of flower morphologies. Our study provides data on how long it takes for B. huntii to forage in open field conditions and will be useful when comparing foraging rates in controlled crop systems.

A homeotic shift late in development drives mimetic color variation in a bumble bee.

Natural phenotypic radiations, with their high diversity and convergence, are well-suited for informing how genomic changes translate to natural phenotypic variation. New genomic tools enable discovery in such traditionally nonmodel systems. Here, we characterize the genomic basis of color pattern variation in bumble bees (Hymenoptera, Apidae, Bombus), a group that has undergone extensive convergence of setal color patterns as a result of Müllerian mimicry. In western North America, multiple species converge on local mimicry patterns through parallel shifts of midabdominal segments from red to black. Using genome-wide association, we establish that a cis-regulatory locus between the abdominal fate-determining Hox genes, abd-A and Abd-B, controls the red-black color switch in a western species, Bombus melanopygus Gene expression analysis reveals distinct shifts in Abd-B aligned with the duration of setal pigmentation at the pupal-adult transition. This results in atypical anterior Abd-B expression, a late developmental homeotic shift. Changing expression of Hox genes can have widespread effects, given their important role across segmental phenotypes; however, the late timing reduces this pleiotropy, making Hox genes suitable targets. Analysis of this locus across mimics and relatives reveals that other species follow independent genetic routes to obtain the same phenotypes.